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Special quantum states exist which are quasiclassical quantizations of regions of phase space that 
are weakly chaotic. In a weakly chaotic region, the orbits are quite regular and remain in the region 
for some time before escaping and manifesting possible chaotic behavior. Such phase space regions 
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We review some methods we have introduced to characterize such regions and find analytic formulas 
for the special states and their energies. 
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I. SPECIAL QUANTUM STATES 

There are many special quantum eigenstates, [more 
precisely classes of eigenstates] which often exist in the 
"quantum chaos" problems that have been considered. 
These states are associated with special regions of phase 
space which are 'weakly chaotic'. By this we mean that 
classical orbits are regular and remain in the region for 
a sufficiently long time, ft may be that almost all the 
orbits are chaotic, but that is a long time property of 
little consequence to the formation of quantum states. 
The regions can be considered as being close to a set of 
nonisolated periodic orbits. These periodic orbits are not 
necessarily orbits of the system under study. They could 
rather be orbits of a 'nearby' integrable reference system. 

This paper reviews some recent work by the authors 
on this topic which is partially published jjj], and 
partially submitted for publication || , . Some further 
examples and details will appear in the thesis of Oleg Za- 
itsev pi . We have introduced a new technique which is in 
some ways more general than those previously used. We 
have also shown how to generalize an older method, the 
Born-Oppenheimer approximation, to new situations. 

Remarkably, only two such kinds of states have been 
clearly recognized in the literature, but there are numer- 
ous others which have been overlooked. This is in spite of 
the fact that the quantum problems giving rise to these 
states have been extensively studied numerically. Even 
in cases where the states dominate the effect under con- 
sideration they have not been mentioned. The reason for 
this oversight is probably as follows: a) the special states 
are high quantum number states and they are rare in the 
sense that they are a small fraction of the ordinary states 
existing in the same energy range, and b) most work has 
concentrated on energy levels and energy level statistics 
rather than wavefunctions. 

The one special state known to most physicists and en- 
gineers, even those not especially concerned with 'quan- 
tum chaos', is Lord Rayleigh's whispering gallery mode 
. Although in most respects this case is very well un- 



derstood, we shall mention a new development here, as 
well as providing a new [or rather old] easy way to get 
the main result. 

There are many phenomena associated with whisper- 
ing gallery modes. The whispering gallery in St. Paul's 
cathedral in London is perhaps the first case to be re- 
marked by a scientist. In 1871, G. B. Airy attempted to 
explain the phenomenon. Lord Rayleigh, in his Theory 
of Sound of 1894 disagreed with Airy and gave, at the 
level of ray analysis, the currently accepted explanation. 
It is amusing that the leading order formula taking the 
wave nature of sound into account is expressed in terms of 
Airy functions. The most recent mention of these modes 
that we have found is by Agam and Altshuler 0. They 
have explained the lack of whispering gallery modes in 
the 'Faraday' experiments of Kudrolli, et. al. || which 
shakes a container of water in the shape of a Bunimovich 
billiard, as due to the fact that significant dissipation 
occurs at the boundary, precisely where the whispering 
gallery modes are localized. 

The second class of states is familiar to most quantum 
chaologists if not to the wider physics community. These 
are the bouncing ball modes. There are some phenomena 
associated with these modes, as well ||. 

The most widespread class of 'bouncing ball' states 
occurs in mixed chaotic systems, that is, systems which 
have stable periodic orbits as well as unstable ones. 
Keller and Rubinow Jlfj| considered, in particular, a 
smooth convex billiard which must in general have a min- 
imum diameter where there is an orbit bouncing perpen- 
dicularly to the boundary at two opposing points. Such 
an orbit is stable and there is a volume of phase space 
near the orbit which is invariant under the Hamiltonian 
flow. It's possible to make a harmonic expansion and 
quantize these states. More generally, there are gener- 
alizations of such 'bouncing ball' states to the neighbor- 
hood of any stable periodic orbit. This requires Ti to be 
small enough that the area of the invariant phase space 
on a surface of section perpendicular to the orbit is of 
order Planck's constant h. 
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Although our technique allows an improvement on the 
harmonic expansion, we shall not consider further such 
states. Since they typically occur in mixed chaotic sys- 
tems they have received relatively little attention from 
the quantum chaos community until recently. 

However, there is a rather spectacular generalization 
of these bouncing ball states to some hard chaotic sys- 
tems not possessing stable orbits. The best known are 
the modes of the Bunimovich stadium corresponding to 
a ball bouncing perpendicularly to the straight sides of 
the billiard. These were discovered by MacDonald and 
Kaufman jTlJ in their pioneering numerical work on that 
system. Similar modes are known in other billiards with 
parallel sides @, @. 

These bouncing ball states have been studied from 
many points of view, including periodic orbit theory [{lif , 
p5[ . Paradoxically, much work has been aimed at 'get- 
ting rid of them'. This is because the classical bouncing 
ball periodic orbits dominate the 'trace formula' and the 
corresponding quantum states dominate the oscillatory 
part of the density of states, thus obscuring the under- 
lying 'quantum chaos' of the energy level correlations, 
which was deemed more interesting. 

Although these special states are associated with peri- 
odic orbits, they are not to be confused with 'scars' |fl6| . 
Scarred states are certainly of great interest. They are 
associated with unstable periodic orbits, and are a rela- 
tively weak effect. The special states are associated with 
a class of reference periodic orbits and are a large effect. 
They are strongly localized in some sense or another, and 
for this reason have the potential to lead to distinct phe- 
nomena. In fact, sometimes one of a sequence of special 
states can be regarded as being scarred. 

The classical limit of the systems under consideration 
may be strongly chaotic, or mixed chaos, or pseudoin- 
tegrable, or even, if you wish, integrable. These are 
long time characterizations of the classical motion, how- 
ever, and have no direct correlation with quantum states. 
To support a special state, the system must be weakly 
chaotic in the sense that there is region of phase space 
where the classical motion is nearly integrable for a short 
but long enough time. The size of this region and the 
characteristic times are functions of Ti. There is also a 
condition on the way that orbits escape from this region. 



II. SOME EXAMPLES 

Here are some examples of systems with special states 
which have not been remarked in the literature. We con- 
sider only two dimensional billiards, for simplicity. We 
begin with three variations on the Bunimovich stadium, 
which has horizontal straight parallel sides of length 2a 
capped by semicircles of radius R, where a/R is of or- 
der unity. We consider the solutions of the Hclmholtz 
equation (V 2 + A; 2 ) ^(x, y) — 0, where vanishes on the 
boundary for values of k 2 on the spectrum. We consider 
only cases where k is large, but not too large. That is, 
we do not automatically take the limit k — ► oo since 



the states of interest will be of negligible number in that 
limit. This is equivalent to considering fi to be small. 



A. Bunimovich stadium 

First, we review the original Bunimovich stadium 
which has a set of nonisolated bouncing ball orbits, 
i.e. those orbits bouncing perpendicularly between the 
straight sides. Such a set of orbits gives an especially 
large contribution to the Gutzwiller trace formula, as is 
well known and we will remind you below. This leads to 
the result that there is a very strong correlation of energy 
levels. In fact, the energy level correlation is dominated 
by the bouncing ball quantum states. These are states 
of the approximate Born-Oppenheimer form, 



where 



$ n (y|x) = sinn7r 



y + R-£ B (x) 
2(R-£ B (X)) • 



(1) 



(2) 



[The Born-Oppenheimer treatment of the Bunimovich 
billiard was first given in reference p7| and is discussed 
further in reference Here £b(x) — 0, \x\ < a, 



£fl(ar) = R - \j R 2 - Qx\ - a) 2 « ~ (\x\ — a) /R for 
|x| > a, (\x\ — a) « R. This makes $(y|x) vanish when 
y is on the billiard boundary. Treating slowly vary- 

ing parameter in <&(y\x), it is seen that tp m (x) satisfies 

- tfjx) + V{x)tl>(x) = E„Mx) (3) 

where V(x) is a sort of square well potential, 

V(x) « n 2 TT 2 £ B (x)/2R 3 (4) 



and the energy 



k 2 - k 2 



{nir/2R'f + E, n 



(5) 



Thus, for large n, tp m (x) is something like sin(m7r(x + 
a)/2a), with m << n. Therefore E m « (mir/2a) 2 . 

The leading term of the energy, (mr/2R) , is a quan- 
tization of the bouncing ball orbits. The next term, is 
a quantization of the transverse motion, motion in the 
potential V(x). If this quantization is carried out semi- 
classically, it too is expressed in terms of the periodic 
orbits of motion in the potential. 

These states have wavenumbers fc„, m ~ n%/2R+0(l). 
This corresponds to the leading term in the trace formula, 
coming from the nonisolated bouncing ball orbits, which 
varies as exp(4ifci?). The length si = 4i? satisfies, ap- 
proximately, ex.p(isi(k n +N,m — k n .rn)) = 1- It is the lead- 
ing element of the length spectrum. The next element, 
s-2 = 8R also satisfies this condition. In the trace for- 
mula it comes from once repeated bouncing ball orbits. 
Thus the bouncing ball states account the most of the 
correlations in the length spectrum which are multiples 
of 4R. 
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We show in Fig. 1 a density plot of ^10,1, as given 
by Eq. ([!]). [Actually, in most cases we show the abso- 
lute square of the theoretical wavefunction.] The wave- 
functions we show are theoretical. Later we make some 
remarks on the accuracy of the theoretical predictions. 




FIG. 1. Bouncing ball mode of the standard Bunimovich 
stadium. Straight side has length 2a, endcap radius is R = a. 
The quantum numbers are n — 10, m = 1. Figs. 1,2,3,4,6,7 
give the stadium shape and a density plot of the square of the 
theoretiacl wavefunction. 

It is, of course, not true that the bouncing ball orbits 
by themselves give the bouncing ball states These 
orbits occupy a set of measure zero in phase space. Nev- 
ertheless, the most natural way to characterize the phase 
space which supports the states is that it is in a certain 
neighborhood of the bouncing ball periodic orbits. 



B. Slightly sloped stadium 

The first variation is the slightly sloped stadium con- 
sidered by Primak and Smilansky In this case the 
end caps have radius i?(l ± e) and the sides, instead of 
being horizontal, have slopes zteR/a. The parameter e 
is small. [To minimize diffraction, we suppose that the 
sides merge smoothly to the end caps with no change of 
slope. The results are given to leading order in e only] 
In this case, the special states are strongly modified even 
if e is quite small. 

In fact, the formulas above hold upon replacing £g — * 
£s where the major new feature is that, for |x| < a, 

£s « -exR/a. (6) 

This gives for V(x) a steep sided well with a sloping bot- 
tom, 

V(x) w -en 2 Tr 2 x/2aR 2 , (7) 

for I a; I < a. The condition that the slope be large enough 
to change substantially the states from the Bunimovich, 
e = 0, case is en 2 / R 2 > 1/a 2 . If this inequality is strong, 
the small m states tp m are approximately Airy functions 
concentrated near x = a. In Fig. 2 we show the state 
*io,2 for e = 0.05, a = R. 




FIG. 2. Special mode of a slightly sloped Bunimovich sta- 
dium. The slope parameter is e = 0.05, and R = a. The 
quantum numbers are n = 10, m = 2. 

One motivation for considering this case is that the 
nonisolated periodic orbits mathematically disappear for 
any finite e no matter how small. However, if e is small, 
there is a region of phase space which remains close to the 
bouncing ball orbits of the original stadium. Thus, there 
are still special states related to the bouncing ball orbits 
of the standard stadium. However, it is complicated and 
not very enlightening to describe these states in terms of 
the periodic orbits of the slightly sloped stadium itself. 

There is a parametrically different dependence on e of 
the states as compared with the energy correlations. We 
see below that if kRe << 1, the leading terms of the trace 
formula, or in other words, the energy correlations due 
to the bouncing ball states, are little modified p5[ . This 
is the natural expectation, since the condition is simply 
that the wavelength is large compared with the change 
of radius of the endcap. On the other hand, we have just 
found that if ka^ft > 1 the eigenstates and energies are 
substantially changed. Thus, if (ka)~ << e << (fci?) -1 , 
the states are changed but the correlations are not. 



C. Baseball stadium 

As e in the previous example becomes of order unity, 
the special state again changes. Note that for finite e, 
there is a set of nonisolated periodic orbits through the 
center of the larger circle, whose angular range is of order 
e. If the endcaps have radii R2, R\ with < R 2 —Ri < 2a, 
the angular range of such periodic orbits is ±9 a measured 
away from the vertical, where s'm9 a = (R2 — R\) /2a. 
This region supports special states ^ n ,m.{r, 9), where r, 9 
are polar coordinates for the large circle. 

The Born-Oppcnheimer approximation used above re- 
lies on having a fast variable and a slow variable. The 
ansatz, Eq. (|l|), was made with y the fast variable and 
x the slow. Fast and slow mean from the classical point 
of view that the motion in the y direction is much faster 
than in the x direction. From a wave viewpoint it means 
that \d$/dy\ >> \d$/dx\ which allows the x dependence 
of $ to be treated parametrically. 

The states of this section would seem to have a fast 
radial motion and a slow angular motion. However, close 
to the center of the circle this separation is not valid. 
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Nevertheless, the Born-Oppenheimer approximate state 
gives the result for the outer region, and that suffices. 

Thus, the states are approximated as above, replacing 
the fast coordinate y — > r and the slow coordinate x — ► 9. 
Then 

$(r|fl) « cos (fcr + K ^ ~ I + a(fl)) (8) 

for fcr large compared with 2/. This function is an 
asymptotic Bessel function of order v, [a superposi- 
tion of Bessel and Neumann functions] and i/j(0) sat- 
isfies f(6) + v{9) 2 iP{9) = 0. The billiard boundary 
is rg(9) = i?2 + £c(0) an( l ^ must vanish for r = 
r a (0). Here £c = 0, -6> a < 6 < tt + 9 a and £ c = 
i? 2 (l/cos(0 + 6> a ) - 1) m ±i? 2 ((9 + 6» Q ) 2 for < -9 a and 
($ + $a) 2 << 1, with a similar formula near 9 = ir + 8 a . It 
can be shown that a = \k{^c{0 + tt) — £c(#)) ± tt/4 and 
^ 2 = £ m - V{0) with V(0) = k 2 R 2 (£, c {6 + ir)+ £ c (0)). 
For large enough k there is an approximate solution 
ip(9) = sm(rmr(6 + a )/26 a ), \9\ < 6 a , ip(9) vanishes out- 
side this region for \6\ < tt, and ip(9) = ±ip(9 + tt). For 
this solution E m = (rrnt/29 a ) 2 . The energy parameter k 
solves kR2 + E m /kR2 = (n=F For large n this has the 

approximate solution k 2 m = ((n qp \)k / R2) 2 — E m jR\. 
There are 2n — | =F \ radial nodes. 

We show this state for n = 10, m = 1, in Fig. 3. Notice 
some changes of sign as compared with the bouncing ball 
between the straight sides. In particular, a deviation of 
the billiard sides in the direction of narrowing the channel 
gives a repulsive effective potential V , thus containing the 
wave function. For a billiard which is a deviation from a 
circle, an outward deviation gives the repulsive effective 
potential. 



studied in this system H]. Take rg = R + £d(0) to de- 
scribe the billiard, with <to(0) small. For the Bunimovich 
stadium with side a << R , £d ~ a |cos0|. Polar coordi- 
nates at the center of the billiard are used and approxi- 
mations neglecting (a/R) 2 have been made. 

In this case the potential V(9) — 2k 2 aR |cos 0\ , an 
attractive triangular well near 8 = ±^7r. It is periodic 
with period tt. For large k 2 aR, the lowest eigenstates will 
be concentrated near 9 = ±\tt. We show the state n = 
10, to = 2 in Fig. 4. 




FIG. 4. Special low angular momentum mode in a Buni- 
movich stadium with a short side, a — 0.17?. Quantum num- 
bers are n — 10, m = 2. 

All of these variations on the Bunimovich stadium 
have been considered in the literature, except for the 
baseball stadium. The only special states previously 
pointed out, however, were for the Bunimovich stadium 
itself. 




FIG. 3. Special mode of a 'baseball' stadium. R2 — 1.5a, 
Ri — 0.5a. Quantum numbers are n — 10, m = 1. The theo- 
retical wave function in this and the next figure is inaccurate 
near the center of the circle. 



D. Almost circular stadium 

The above Eq. (||) can be applied to find the low angu- 
lar momentum states of any almost circular billiard. En- 
ergy levels and wavefunction statistics were extensively 



III. FINDING AND APPROXIMATING SPECIAL 
STATES 

A. History 

The whispering gallery modes and the bouncing ball 
modes corresponding to stable periodic orbits were first 
quantized by Keller and Rubinow fic|| by what is some- 
times called the ray method. This technique starts with 
the classical mechanics or rays, and exploits caustics and 
adiabatic invariants. However, the assumptions usually 
made are unnecessarily strong, and these special states 
still exist even if the caustics are only a short time ap- 
proximation. 

Another group of methods starts from the partial dif- 
ferential equation, the Helmholtz equation in the case of 
billiards. The parabolic e qua tion method invented inde- 
pendently by Leontovich flS|| and by Fock []l9| chooses 
coordinates astutely, and finds appropriate scale fac- 
tors, allowing an approximation to the partial differen- 
tial equation. It often relies on the ray method to moti- 
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vate the manipulation of the PDE. The etalon method of 
Babich and Buldyrev [ ^p| is an improvement of this which 
involves choosing a characteristic example, or 'etalon', 
which captures the essential features of the given prob- 
lem. The extended Born-Oppenheimer method [EBO] has 
been used in this context only for bouncing ball states 
between parallel sides of a billiard. We showed above a 
couple of generalizations of this technique. Other gener- 
alizations are given elsewhere || . 

The last three methods are closely related. They use 
different language and motivation and make approxima- 
tions in a different order, but the main result is the same. 
Systematic corrections to the leading result have been ex- 
tensively studied for certain examples, and these correc- 
tions are apparently of a different form for the different 
methods. We prefer the EBO when it applies since it is 
simpler and better known than the other techniques. 



w here L is the distance between boundary points, k = 
\/~2mE/h~ and /i = 7r for Dirichlet conditions. We shall 
suppress /j, and write the prefactor as (.) for simplicity. 
The essential variation of T is given by the exponen- 
tial dependence on kL = HkL/h = S(s, s'\E)/h~, where S 
is the action of the classical straight line orbit between 
boundary points s, s' . 

Because h is supposed to be small compared with typ- 
ical classical actions the exponential is rapidly varying, 
and the method of stationary phase usually applies to in- 
tegrals in which T appears. We deal always with billiard 
examples for which it is convenient to take fi = 2m = 1, 
and kL is typically large. 

The action S generates the classical surface of section 
map from 

Ps = dS/ds, Ps , = -dS/ds'. (11) 



B. Bogomolny operator 

We also have introduced a technique based on Bogo- 
molny's surface of section transfer operator |^lf which in 
certain ways is more general than the above methods. 
It also has the advantage that the transfer operator T is 
closely related to the trace formulas and to important 're- 
summations' of the trace formulas. In addition, it makes 
somewhat more precise the notion of a region of phase 
space which is nearly integrable for short times. 

The operator T(s, s'\E) introduced by Bogomolny is a 
generalization of the boundary integral method used to ob- 
tain numerical solutions for billiard problems. The main 
equation of the boundary integral method is derived us- 
ing Green's theorem. It takes the form 



M (s)= / ds'K( S , S '\E)n(s'), 
Jd 



(9) 



a Fredholm integral equation. The kernel or operator K 
depends parametrically on the energy E. Eq. (^) has 
a nontrivial solution only when E is on the spectrum. 
For Dirichlet conditions, /i(s) = d^(r)/dn, the normal 
derivative of the wavefunction evaluated at the position 
on the boundary labelled by s the distance along the 
perimeter. The integral in Eq. ([)]) is over this boundary. 

The boundary coordinate s together with its conju- 
gate momentum, the momentum parallel to the bound- 
ary p s , gives the Birkhoff coordinates for a surface of 
section of a billiard. Bogomolny's method essentially ap- 
proximates K(s,s'\E) by its asymptotic form, which we 
call T(s, s'\E). It generalizes to a broad class of surfaces 
of section, not restricted to billiards, for which there are 
often also exact kernels. For a billiard, in Birkhoff coor- 
dinates, 



T(s,s'\E) = 



' k\d 2 L(s,s')/dsds'\ \ 2 j(fci(s v) +A1 ) 



2m 



Composition of powers of T have intermediate points of 
stationary phase determined by this map so in effect, 
longer and longer orbits can be built up by iteration of 
the T operator. 

Corresponding to Eq. (Jot) is the fundamental equation 



tp(s) = / ds'T(s,s'\E)^(s'). 



(12) 



Our new method, for the cases corresponding to special 
states, solves this equation directly and quasiclassically, 
that is, it exploits the stationary phase approximation to 
do the integral. 

Eq. (|l^) has solutions only for E values on the quasi- 
classical spectrum such that 



D(E) = det(l - T(E)) = 



(13) 



This Fredholm determinant is a well defined version of 
the dynamical zeta function [p2| , a resummation of the 
trace formula p3|] . The trace formula itself is given by 



, . -1 dlnD(E) 

d osc (E) = — Im — 

7T dE 



-1 dTrln(l-T(E)) 
= Im 

TT dE 

TT dE ^ r 



(14) 



(10) 



Here the density of states, d(E), is expressed as 

d(E) = ^S(E- E a ) = d Weyl (E) + d osc (E) (15) 



where dweyi is the smoothed density of states. The traces 
of powers of T are expressed in terms of periodic orbits 
when the integrals are done in stationary phase. In the 
case of the special states, as well as the integrable case, 
however, only r — 1 integrals are well approximated by 
that method. 
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C. Special states and periodic orbits 



D. Nearly integrable phase space region 



We take as an example the slightly sloped billiard de- 
scribed above. Introduce a reference billiard system con- 
sisting of a channel of width 2R, i.e. with boundaries 
at y = ±R. The actual orbits are compared with the 
bouncing ball orbits of this channel. 

We want to choose a surface of section such that the 
nominal periodic orbits are described by a diagonal ele- 
ment, T(x, x). Such a surface of section is the upper half 
of the billiard. It is convenient to label position along it 
by the corresponding position x on the upper side of the 
reference billiard. [Note that a simple change of coordi- 
nates, e.g. s — > s(x) does not affect the result.] 

We now approximate the actual T(x, x') in a way valid 
for | a; — x'\ « R. Thus, we may expand 

L(x, x') « AR + ^ (x - x'f - ts(x) - tstf) (16) 

where £<?, given by Eq. ([|) above, is small except in 
the endcaps. The first two terms in this formula are an 
approximation to the perfect channel, and £s approxi- 
mates the difference between the sloped billiard and the 
channel. 

Assuming for the moment that k!;s( x ) << 1, we can 
approximate 



ikR 



(4+i(^) : 



(17) 



T(x,x')^(.)e 

x e -*k(£s(*)+6s(a:')) 
raQe i*«(4 + i(*#) 2 ) 

x(l-tt(&(a;) +&(*')))■ (18) 

This T operator has a form that allows direct solution 
of Eq. (T^), Next assume that the solution tp{x) of 
Eq. ( |l2| ) is slowly varying and expand ip(x') rs + 
(x' — x) ip'(x)+^ (x' — x) 2 ip"(x). The rapidly varying ex- 
ponential in T ensures that the important contributions 
to the x' integral lie close to x, according to the estimate 
(x'-x) 2 ~R/k. 

One may carry out the integral of Eq. ( |l2|) and obtain 
the equation tp(x) — exp(i(4fci? — E m R/2k))ip{x) with 
the condition that ip satisfies the Schrodinger equation 
(||) with potential V of Eq. (Q), where the replacement 
nir — > 2k R is made. The energy eigenvalues k^ rn are 
found by insisting that the exponential is unity. Thus, 



^k nm R E m R/2k nm — 27T7Z 



(19) 



Since fc£g becomes large as x goes into the endcap re- 
gion, the approximation of Eq. (|18| ) seems to break down. 
However, the wave function becomes small in that region. 
As long as the approximate wave function continues to 
be small, the approximation is satisfactory. 

It is also possible to extend this result to a parameter 
regime such that fc£s becomes of order unity or larger |Q . 
It is only necessary that fc£s be slowly varying compared 
with the leading term. 



The above formulation can be regarded as defining a 
phase space region where the system is nearly integrable. 
The two arguments of T together with the surface of sec- 
tion map give an area on the 2-dimensional phase space 
such that T is well approximated by Eq. ([l8|). This, in 
turn, is the T operator for an integrable system, char- 
acterized by being a function of the coordinate differ- 
ence, together with a leading correction. For the case 
discussed, the surface of section phase sp ace i s approxi- 
mately |x| < a, \p\ < hk \x — x'\ /R w hy/k/R. The two- 
dimensional surface of section phase space region defines 
a four dimensional phase space region, from which the 
orbits only slowly escape. 

A classical orbit at a typical point in this phase space 
region is not too different from a nearby classical orbit 
of the ideal integrable region. Of course, continuing this 
orbit for long times allows it to escape from the region. 
In the case at hand, almost all such orbits are chaotic 
with positive Lyapunov exponents. 

However, if the fundamental equation Tip — ip has a 
solution ip(x) which is large only in the nearly integrable 
region, there is no need to consider the very long orbits. 
That gives another condition on the perturbation. In 
terms of the effective potential, V(x), it means that V 
must become repulsive outside the region. The sign of V 
depends on the sign of the perturbation, and on the sign 
of the quadratic term in the integrable part of T. Thus, 
for bouncing ball states between parallel sides, a per- 
turbation shortening the bounce path is repulsive, while 
for radial bounces between circular sides, a perturbation 
lengthening the path is repulsive. 



E. Trace formula 



Consider the contribution to the trace formula from 
the period 1 orbits, of Eq. ©, d oscA = ^Im^TrT. In 
our case, TrT = JdxT(x, x). The hard chaos case usu- 
ally considered assumes that this integral can be done 
by stationary phase, with contributions only from the 
neighborhood of one or more points. Using Eq. jl^ ) we 
see that the stationary phase approximation is not good 
unless eka » 1. In the case of the opposite strong in- 
equality, eka << 1 the x integral gives just a factor of 
2a, to leading approximation. However, we saw that if 
yfeka > 1, there are strong effects on the wave functions. 
Thus, although the wave functions are sensitive to the 
slope, the energy level correlations are not in this param- 
eter regime. 

This is a fairly general result. Namely, special states 
are quite sensitive to perturbations. The energy level 
correlations are much less sensitive to the same pertur- 
bation. 
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IV. MORE EXAMPLES 

We briefly give a few more examples. The problems al- 
ready mentioned can be solved by the Born-Oppcnhcimcr 
method. In the examples of this section, the Born- 
Oppenheimer method is more limited or more cumber- 
some, or impossible to use, but the Bogomolny operator 
method succeeds. 

First, the whispering gallery modes can easily be stud- 
ied. These modes are concentrated near the boundary of 
a sufficiently smooth and sufficiently convex billiard. The 
case in which the corresponding classical orbits are also 
confined near the boundary is the only one studied in de- 
tail in the literature. Assuming a smooth enough convex 
billiard with everywhere positive curvature on the bound- 
ary, the appropriate coordinates are s, the distance along 
the perimeter, and p, the distance from the perimeter to- 
ward the center of curvature at s. Let R(s) be the radius 
of curvature at point s. The Born-Oppcnhcimcr ansatz is 
*(r|«) = e* ks <5>(p\ S )iP(s), Hp\s) = C(s)J v{s) (k(R(s)-p)) 
and ip(s) = e~ tkR ( s >f( s K Here J v is a Bessel function of 
the first kind, with variable index v and C(s) is a slowly 
varying prefactor, determined by normalization. 

The whispering gallery modes are characterized by 
large v : slightly less than kR. It turns out that v — 
kR(l — /) where / is small. One choice is using the 
Born-Oppenheimer approach, treating the s variation as 
slow compared to p, [except for the explicit factor e lks ], 
and determining v(s) by the condition J v t s \(kR(s)) = 0. 
[A standard asymptotic formula expresses J„ in terms of 
an Airy function.] Alternatively, if) can be determined as 
the solution of Tip = ip and $ can be obtained from that. 
The result is f(s) — z m /2s (fc_R(s))3, where z m is a zero 
of the Airy function. 

It turns out that the condition \d$(p\s)/dp\ >> 
\d&(p\s)/ds\ that the Born-Oppcnheimer approximation 
is valid is the same as the classical condition for a caus- 
tic. [The quantizing caustic would be at approximately 
p = R(s)f(s).] The method using the T operator can 
be generalized, however, to study cases where there are 
points of vanishing curvature on the boundary and caus- 
tics do not exist. Associated with a caustic is an adiabatic 
invariant. As an orbit passes a zero curvature point, the 
adiabatic invariant jumps to a new value. The T opera- 
tor can account for this, but the Born-Oppcnheimer wave 
function cannot be used. We show in Fig. 5 a whispering 
gallery state in a stadium billiard with small short side 
a = Q.05R. There are no caustics in this case. 

A second example is the slightly sloped trapezoidal 
billiard of Morse and Feshbach [53 , which was recently 
reconsidered by Kaplan and Heller |^6| . The remarkable 
states shown in the next two figures were not mentioned, 
however. As shown in Fig. 6, this billiard is almost a 
square, with vertical sides, say, x = 0, x = 1, and hor- 
izontal sides y = 1, y = ex, where < e << 1. The 
Born-Oppenheimer method may be used to find states 
concentrated near x — 0, similar to the states of the 
slightly sloped stadium. However, there are states asso- 
ciated with the (1,1) period orbits of the unperturbed 



square. These orbits are rectangles making a 45° angle 
with the a;— axis. This may be solved using the Bogo- 
molny technique by extending the billiard antiperiodi- 
cally in the x direction, and using as surface of section 
the upper boundary, y = 1. The effective potential is 
V(x) cx ek 2 \x\ , —1 < x < 1, repeated with period 2. 
States of this extended problem are superimposed to find 
a solution. A typical case is shown in Fig. 7. 




FIG. 5. Contour plot of a numerically obtained whispering 
gallery mode in a stadium billiard, a = 0.05-R. There are no 
caustics in this case, and almost all classical orbits circulating 
near the boundary eventually escape. The distance between 
the short parallel lines is the length of the straight side. The 
wavenumber is k — 242.761 1/R. 

A last example is again the almost circular stadium. 
We now study states near a higher period orbit. In Fig. 
8 we show a numerical state near the (1,4) almost square 
orbits in a stadium with side a/R = 0.01. We have not 
found any convenient coordinates which we can classify 
as fast and slow, so the EBO does not work. The T— 
operator approach is straightforward, but a little compli- 
cated [Q. In the T— operator approximation, there are 
two nearly degenerate states. That is, the T operator has 
inequivalent solutions of the same energy corresponding 
to orbits moving clockwise or counterclockwise. These 
states do not couple in the T operator approximation, 
and the appropriate symmetric combinations have nearly 
the same energy. 

From this figure it is seen that the states are made 
up of waves along the rays of classical angular momen- 
tum I w hk/V2, i.e. along straight lines whose closest 
approach to the center is l/v2- A perfect circle would 
have a caustic of this radius. Because the stadium is 
chaotic, such a caustic does not exist, except as a short 
time approximation. 

V. SUMMARY 

We have reviewed some of the salient features of spe- 
cial classes of states which often occur in the systems 
of interest to quantum chaologists. Perhaps these cases 
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are common because there is a tendency to construct bil- 
liards of straight lines and circles. The system as a whole 
may be chaotic, in other words, almost all orbits of the 
system may have positive Lyapunov exponents. How- 
ever, for some shorter time, a subset of phase space may 
be close to that of a set of nonisolated periodic orbits of 
some integrable reference system, and this leads to the 
special states. 




FIG. 6. State along the long side of a nearly square trape- 
zoid. Slope of bottom is e = 1/16. Quantum numbers are 
n = 22, m — 2. In this figure and the next, the dashed line is 
the a;— axis. 

These special states have rather striking wavefunc- 
tions. As a result, there are phenomena and even possi- 
ble applications associated with them. In particular, the 
whispering gallery modes have long been known to give 
rise to interesting effects. One standard 'application' is 
that special states often appear in weakly perturbed in- 
tegrable systems. If these states are to be avoided, a 
regular resonant cavity must be constructed much more 
precisely than the condition 5x « A, where 5x is the 
deviation from the ideal. Also, because there are a se- 
quence of special states regularly spaced in energy, the 
special states often numerically dominate the trace for- 
mula. 

Typically, there are two or more scales of variation in 
connection with the special states which can be identi- 
fied. If coordinates can be found such that one coordinate 
is fast and the other is slow, standard adiabatic approxi- 
mations, such as the Born-Oppenhcimcr approximation, 
can be used. 

Bogomolny introduced a surface of section transfer op- 
erator some time ago as a means of studying the spec- 
trum, i.e. the trace formula. We have shown how this op- 
erator can also used to find the special eigenstates. The 
technique works if the operator can be approximated as a 
rapidly varying integrable part, and a more slowly vary- 
ing correction to it. It is in some ways more general than 
the other methods. 




FIG. 8. A contour plot of a numerically obtained state near 
the (1,4) periodic orbits in a Bunimovich stadium billiard with 
a short straight side, a/R = 0.01. The wavelength of the state 
is A = 2.23969a. The parallel lines show the length of the short 
side, 2a. 

The special states are typically rather rare, in the sense 
that they are a small fraction of all the states in a given 
(high) energy range. In leading approximation, they do 
not couple to the other states. The energies of the states 
are predicted to good approximation, relative to the en- 
ergy spacing of the given class of states, and even better 
absolutely. However, the accuracy is not necessarily good 
compared with the mean level spacing of all the levels. 
Further, it may happen that 'accidentally' there is a non- 
special state with energy very close to that of the special 
state. Then terms neglected in our approximation can 
mix these states. Many phenomena are independent of 
such mixing, however. 
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In this article we gave a number of examples of such 
special states, several of which appear for the first time in 
print. We hope our pictures will tempt an experimental- 
ist to find some of these states in one of the several sys- 
tems, water trays, acoustic plates, microwave and laser 
cavities, optical fibers, quantum dots, . . . , to which the 
theory applies. 
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